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Polaroids:  a  new  cool  in  non- convex 

and  in  integer  programming 


by 

Claude-Alain  Burdet 

ABSTRACT 


This  paper  presents  a  generalization, called  Polaroid ^of  the  concept 
of  polar  sets. 

A  list  of  properties  satisfied  by  polaroids  is  established  indicat¬ 
ing  that  the  new  concept  may  be  fruitfully  used  in  an  area  of  non-convei: 
(called  here  polar)  programming  as  well  as  in  integer  programming,  by 
means  of  Polaroid  c.its;  this  class  of  new  cuts  contains  the  ones  defined  by 

Tuy  for  concave  programming  (a  special  case  of  polar  programming)  and  by 
'  % 

Balas  for  integer  programming;  it  furthermore  provides  for  new  degrees 
of  freedom  in  the  construction  of  algorithms  in  the  above-mentioned  areas 
of  mathematical  programming. 


1. 


1)  Preliminaries 

In  this  introductory  section  we  present  the  definitions  and  some 
relevant  properties  of  a  new  concept:  Polaroid  sets  and  functions; 
these  mathematical  objects  are  derived  directly  from  the  theory  of  con¬ 
vex  sets  and  represent  a  generalisation  of  polar  sets;  historically, 
polar  relationships  have  been  one  of  the  main  topics  in  projective  geome¬ 
try  where  the  involutory  correspondence  of  poles  and  polars  is  of  prime 
interest.  The  first  generalization  to  polar  sets  can  be  found  in  Min¬ 
kowski  [8  J;  these  sets  have  since  played  an  increasingly  important 
role  in  the  area  of  convex  analysis  and  mathematical  programming  (see 
for  instance,  the  treatise  [10]  by  Rockafellar) ;  more  recently,  Salas 
[2]  has  uncovered  an  interesting  application  of  polar  sets  to  integer 
programme ,  there  the  use  of  a  positive-definite  quadratic  form  (n- 
dimensional  sphere  or  ellipsoid)  allows  him  to  define  his  outer-polar 
sets  which  enjoy  all  the  desirable  properties  for  the  convex  outer- 
domain  theory  of  valid  intersection  cuts  (further  aspects  of  this  theory 
have  been  discussed  by  Glover  [5]  and  Kurdet  [3,4]). 

Our  objective  here  is  merely  to  present  some  fundamental  properties 
of  polaroids,  and  to  indicate  how  they  can  be  fruitful  in  an  area  of 
non-convex  linearly  constrained  programming  (called  here  polar  program¬ 
ming)  as  well  as  in  integer  (primarily  zero-one)  programming. 

When  the  constrained  set  is  polyhedral,  the  use  of  polaroids  can  be 
viewed  as  a  generalization  of  the  intersection  approach  of  Hoang  Tuy  [7] 
for  concave  programming  or  of  the  intersection  cut  approach  initiated  by 
Balas  [  1  ]  in  integer  programming.  The  generalization  with  respect  to 
[ 7  ]  is  that  Polaroid  cuts  can  be  defined  for  non-concave  problems; 
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furthermore,  for  concave  problems  Tuy's  cuts  are  uniformly  dominated  by 
Polaroid  cuts;  indeed  Tuy's  approach  can  be  imbedded  as  a  trivial  special  case 
in  the  present  r.heory.  With  respect  to  the  results  of  Balas  [3,2],  the 
generalization  consists  in  the  fact  that: 

a)  convex  polaxoids  can  be  defined  for  a  more  general  class 
of  functions  than  positive  definite  quadratic  forms  as 
in  [2  ];  adequately  constructed  polaroids 

can  be  used  as  convex  outer-domains  to  generate  new 
cutting  planes. 

b)  the  convexity  requirement  imposed  on  outer-domains 
merely  plays  the  role  of  a  sufficiency  condition  in  the 
construction  of  an  intersection  cut^ 

and  it  can  be  relaxed;  it  is  shown  that  non- 

convex  polaroids  may  very  well  be  used  to  generate  valid 

Polaroid  cuts. 

Men  culpa:  Tne  underlying  goal  of  this  report  is  not  primarily  focused 
on  numerical  and  computational  aspects  of  the  uncovered 
properties;  the  basic  idea  was  to  start  with  as  general  an  object  as  pos¬ 
sible  (viz.  Polaroid  sets)  and  to  test  those  properties  which  seem  prom¬ 
ising  for  the  global  optimization  of  non-convex  problems.  The  analysis 
resulted  in  a  hierarchy  of  properties  ranked  by  increasingly  strong 
assumptions.  In  a  second  effort  it  was  found,  however,  that  many  of  the 
additional  assumptions  leading  to  the  more  sophisticated  (higher  ranked) 
properties  (such  as  convexity)  certainly  are  convenient  because  they  pro¬ 
vide  for  automatic  sufficiency  conditions,  but  at  the  same  time  they  seem 
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to  be  unnecessarily  tamin6  the  power  of  the  approach,  ^or  instance,  for 
the  construction  of  a  cutting  plane  it  is  believed  (as  indicaied  in  the  con¬ 
clusions)  that  a  fruitful  line  of  research  (fn  order  to  attain  computational 
efficiency)  would  be  to  proceed  along  the  following  way: 


By  convention  set  P*(k)  =  0  ,  Vk  <  min  {f(x,y)  {  xeP] 

y 

Theorem  1  (Inclusion  theorem): 

For  any  closed  sets  P  and  Q  C]Rn  ,  one  has  the  following  implica¬ 
tions: 

If  Q  c  P  then  Q*(k)  z>  P*(k)  ;  and 


If  kx  >  k2  then  P*(k^)  =5  P*(k2)  . 


4. 


Proof:  The  assertions  follow  from  (1): 

Q*(k)  =  {y I  f (x,y)  <  k  ,  VxcQ)  3  (y |  f(x,y)  <  k  ,  VxcP  =>  Q}  *=  P*(k) 

P*^,)®  { y!  f(x,y)  <  k1#  VxePj  (yj  f(x,y>  <  k2  <  ^  ,  Vxep}  =  P*(k2> 

Q.E.D. 

Theorem  2  (Union  theorem):  For  any  closed  sets  P  ,  Q  one  has 

(P  U  Q)*  (k)  «  P*(k)  0  Q*(k) 

Proof:  (P  U  Q)*  (k)  -  (y|  f(x,y)  <  k  ,  Vx«(P  U  Q)}  • 

»  (y|  i'(x,y)  <  k  ,  Vx«P}  n  {y|  f(x,y)  <  k  ,  VxcQ} 

“  P*(k)  f  Q*(k)  .  Q.E.D. 

Corollary  2.1:  Let 

(P  A  Q)  «  cl  {x  |  x«(P  U  Q)  *nd  xrf(P  0  Q)} 

Then  one  has 

(p  n  q)*  (k)  n  (P  a  Q)*  (k)  -  p*(k)  n  g*(k) 

Proof:  (P  U  Q)*  (k)  -  ( (P  U  Q)  U  (P  A  Q)]*(k)  -  (P  0  Q)*(k>  0  (P  A  QMk) 

On  the  other  hand,  theorem  2  implies: 

(P  U  Q)*(k)  »  P*(k)  n  Q*(k)  Q.E.D. 

Theorem  3  (Intersection  theorem): 

(P  n  Q)*  (k)  3  P*(k)  U  Q*(k) 

Proof:  (P  n  Q)*  (k)  =  {y  |  f(x,y)  <  k  ,  Vxe(P  D  Q)} 

3  {y  I  f (x,y)  <  k  ,  VxeP}  U  {y  J  f(x,y)  <  k  ,  VxeQ]  = 
-  P*(k)  U  Q*(k)  . 


Q.E.D. 
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Corollary  2.2: 


P*  (k) 
Q*(k) 


(p  n  q>*  (k)  n  (p  i  (p  r  q»*  (k> 
(p  n  q>*  (k)  n  (Q  a  (p  n  q>>*  (k) 


Proof:  P  *  (P  D  Q)  U  (PA  (P  f!  Q))  ,  and  similarly  for  Q  ;  applying 

the  union  theorem  completes  the  proof.  Q.E.D. 

Corollary  2.3: 

P* (k)  u  Q*(k)  =  (P  n  Q)*(k)  n  [(Pi  (P  n  Q))*(k)  u  (Q  A  (P  n  Q))*(k)J  . 

^  (P  n  Q)*(k)  n  (PA  Q)*(k)  =  P*(k)  n  Q*(k) 

Proof:  The  (first)  equality  is  obtained  directly  from  the  corollary  2.2; 
the  inclusion  can  be  dej ived  by  inspection  of  the  polaroid  sets 

"ft  "*C  "jc 

in  the  bracket  or  directly  from  corollary  2.1,  since  P  UQ  ^P  HQ  . 
Definition  2:  The  function  g  =  g(x)  =  f(x,x)  is  said  polarized  by  f  . 

Denote  the  level  set  of  g  by  lev^  g  =  {x|  g(x)  <  a}  . 


Defi-ltion  3:  An  arbitrary  (closed  compact.)  set  S  satisfying 
0  *  S  ft  P  c.  ievk  g 

is  called  a  valid  cut  at  the  level  k  . 


Theorem  4:  A  cut  S  f  valid  at  the  level  f  is  valid  at  all  higher 
levels  k2(>  kj)  . 

Proof:  Immediate  since  lev  g  <=  lev  g  by  definition. 

1  z 


Q.E.D. 


givaa  nz-v  w 
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3)  Convex  Polaroids 

Theorem  5:  The  polaroid  P*(k)  is  convex  Vk 

iff  f(x,y)  is  guasi-convex  in  y,  for  all  xgP  . 


1  2 

Proof:  Let  y  and  y  be  arbitrary  points  in  P*(*)  ;  and  set 

y3  =  Xy1  +  (l-X)y2  ,  0  <  X  <  1 
then  VxeP  ,  one  has 

f(x,y3)  =  f^Xy1  +  (l-X)y2)  <  max  (fCx.y^)  ,  f(x,y2)}  =  k 

Furthermore,  assume  P*(k)  convex  for  all  k  ;  thus  let 

3  12  3 

y  =  Xy  +  (l-X)y  ,  then  one  has  y  e  P*(k)  ,  for  all  k  such  that 

1  0 

y  and  y  e  P*(k)  i.e.  for  all  k  such  that: 
k  >  max  {f(x,y  )  ,  f(x,y  )  j  xeP|  • 

In  particular  for 

k  =  k  max  (f(x,y^)  ,  f(x,y2)  j  xep)  . 

3  it  • 

Since  y  e  P  (k)  one  has  Vxep  : 


f(x,y3)  <  k  =  max  {f(x,y1)  ,  f(x,y2)}  ,  VxeP 

Q.E.D. 

In  order  to  acquire  at  this  point  a  better  geometrical  feeling  for  the  con¬ 
tent  of  the  statement  in  theorem  5  let  us  review  some  classical  results: 

n 

1)  Let  f  be  the  euclidean  scalar  product  f(x,y)  =  £  x.y.  =  <x,y> 

i=l  1  x 

which  yields  (the  square  of)  the  euclidean  norm  as  polarized 

n  2 

function  e(x)  =  £  x  . 

i=l 
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Following  the  classical  definition  [10]  one  finds 
that,  in  this  case,  the  polaroid  P*(l)  is  nothing  but  the 
polar  set  P* 

P*(l)  =  P*  =  {x*  {  <x,x*>  <  1  ,  VxeP] 


In  [2]  Balas  introduced  a  generalization  of  the  classical 
concept  by  allowing  for  a  scaling  factor  k  in  the 
above  definition.  Clearly  this  can  be  absorbed  in  our 
definition,  either  by  considering  the  polaroid  P  (k) 
with  respect  to  the  same  polaroid  function  f  ,  or  by 
changing  the  polaroid  function  to  f  =  —  f  and  retaining 
the  parameter  1  :  P  (1)  . 

Since  f  is  bilinear  it  is  quasi-convex  in  y  and  theorem  5 
paraphrases,  in  this  case,  the  convex  property  of 
polar  sets. 

2)  Let  £  be  the  scalar  product  function  corresponding  to  a 
general  Riemann  metric,  with  arbitrary  (real  symmetric) 
metric  tensor  ;  i.e. 


n 


£<X,y)  =  ■  L,  8ikVk  ’ 

i  ,  k— 1 

for  the  polarized  function  g  one  has 


n 


8(10  =  ■  L  gikVk 

i,k=l 


Since  the  parameter  k  can  always  be  absorbed  by  the  tensor 
g.^  ,  the  polaroid  P*  =  P*(i)  is  a  genuine  generalization 
of  the  polar  concept,  to  arbitrary  quadratic  forms  (not 
necessarily  definite  or  positive). 

Since  f  is  again  bilinear,  one  finds  from  theorem  5 
that  the  Polaroid  P*  generated  from  an  arbitrary  set  P 
is  convex;  in  a  follow-up  paper,  we  show  how  the  Polaroid 
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theory  can  be  used  to  solve  an  arbitrary  non-convex  quad¬ 
ratic  programming  problem  (particularly  in  the  indefinite  case)  [1?]# 


3)  The  following  extension  of  a  Riemann  metric  was  found 
useful  to  generate  convex  outer-domains  in  integer  pro¬ 
gramming  [h  ] . 

n 

Let  f(x,y)  =  E  g.x  y 
i=l  x 


»  /  +.-2 
"  W 


where  (  =0 


for  x^  >  0 


for  a.  =  0  or  aJ  =  0 
x  i 


=  VaI>'2 


for  x.  <  0 
x 


n  +  - 

with  E  A.  =  1  ,  A.  >  0  ,  (a.  and  a.  being  given  quantities), 

x  l  —  i  i 

i=i 


The  corresponding  polarized  function  g 


n 


g(x)  =  E  g  x~  is 
i=l 


piecewise  quadratic,  i.e.  it  is  quadratic  in  each  of  the 
2°  "orthants  of  Fn."  But  here  agaii  because  f(x,y)  is 
linear  in  y  for  every  x  in  a  given  orthant,  the  pola- 

JU 

roid  set  P  (k)  can  be  shown  to  be  convex. 


4)  Consider  a  quasi-convex  function  g(x)  and  define 
f(x,y)  =  g(y)  ,  Vx 

The  polaroid  P*(k)  then  merely  reduces  to  P*(k)  =  lev^  g 
and  theorem  5  restates  the  known  convexity  property  of 
leVj_  g  .  In  general,  however,  there  belongs  many  other 
possibilities  for  choosing  f  =  f(x,y)  such  that  f(x,x)  =  g(x)  ; 
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hence  there  will  correspond  other  polaroid  sets  to 
the  same  polarized  function  g(x)  ;  example  1  illus¬ 
trates  this  situation,  for  instance. 

5)  This  list  could  extend  indefinitely!  As  a  final  example 


consider 


f  =  f(x,y)  = 


T. 

i»j=l 


(cij  +  iAftV  Vj 


and  suppose  that  P  has  the  property  that  the  matrix 
[C  +  Dx]  is  positive  semi-definite  VxeP  .  Then 
f(x,y)  is  convex  VxeP  and  the  polaroid  P*  is 


convex. 


Until  now  we  always  considered  k  as  an  accessory  parameter  whose 
value  played  no  essential  role.  We  now  analyze  polaroids  (convex  or  not) 
corresponding  to  particular  values  of  k  . 

Definition  4:  The  polaroid  P*(k)  defined  by  (1)  is  called  complete  if 

P  p*(k) 

Definition  5:  A  function  f  =  f(x,y)  is  called  complete  on  P  at  fee  level  k  if 

f(x,y)  <  k  ,  Vx,yeP 
Theorem  6  (Coinpletei  css  theorem) : 

The  polaroid  P*(k)  is  complete  iff  f  is  complete  on  p  (at  the  level  k)  . 
Proof:  Taka  any  x'eP  ;  then  f(x,x')  <  k  ,  VxeP  by  hypothesis> when 
f  is  complete;  but  one  has  P*(k)  =  {y}(f(x,y))  <  k  ,  VxeP]  , 


which  shows  that  x'  must  belon  ;  to  P*(k)  and  hence  P  c  P*(k)  . 
Conversely,  suppose  P*(k)  complete;  the  argument  is  by  contradiction: 
suppose  there  exists  a  pair  x,yeP  such  that  f(x,y)  k  ;  in  this 
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case  one  has  either  x^P*(k)  or  y^P*(k)  ,  or  both;  in  any  case 
P  ?  P*(k)  which  is  contrary  to  hypothesis. 

Q.E  .0 . 


The  following  theorems  will  be  useful  in  establishing 
optimality  conditions  for  polar  programming  as  well  as  for  testing  the 
validity  of  polaroid  cuts  in  integer  programming. 

k 

Corollary  6.1:  If  P  (k)  is  complete  then  P  c  lev^g  • 

Proof:  By  hypothesis  one  has  f(x,y)  <  k  ,  Vx,yeP  ;  hence,  in  par¬ 
ticular  for  y  =  x  :  f(x,x)  =  g(x)  <  k  . 


Q.E.D. 


Theorem  7 :  If  f  is  symmetric,  i.e.  f(x,y)  =  f(y,x)  then 

P  c  (P*(k)),r(k) 

Proof:  P  (k)  =  {y|f(x,y)  <  k  ,  Vx^Pl  ,  hence 

f(x,y)  <  k  »  VxgP  ,  Vye?'(k);  but  by  hypothesis 
f(x,y)  =  f (y,x)  <  k  ,  hence  VxeP  ,  one  has 
xe(P  (k))  (k)  =  {y|f(y,u)  <k  ,  VyeP*(k)} 

and  therefore  P  C  (?"(k)r(k)  Q.E.D. 

k 

Corollary  7,1:  If  P  (k)  is  complete  (and  f  symmetric)  then 

p  c  (P*(k))*(k)  C  P*(k) 

k 

Proof:  pep  (k)  implies 

(p'(k))~(k)  =  {y|  f  (x,y)  <  k,  Vx sP*  (k) )  C 

c  {y| f (x,y)  <  k,  Vx€P}  =  P*(k) 


Q.E.D. 


11. 


Let  us  now  consider  the  "boundary"  sets 


bd  P*(k)  =  {yeP*(k)  |  f(x,y)  =  k  ,  for  some  xeP} 
bd  (lev^g)  =  {x  }  g(x)  =  k}  . 

Theorem  8:  Suppose  P*(k)  is  complete,  then  one  has 


[P  0  bd(lev  g)]  c  [p  n  bd  ?*(k>] 

4C 

Proof:  First  let  us  note  that  the  assertion  is  trivial  for 
[P  0  fcd(levkg)]  =  0  . 

Take  xe  [P  D  LdClev^g)]  ,  i.e.  xeP  with  g(x)  =  k  .  Since 
P*(k)  is  complete,  one  also  has  xe?*(k) .  Ihus  xeP*(k)  and 
g(x)  =  f(x,x)  =  k  ,  with  xeP  which  completes  the  proof. 

Q.E.D. 

Theorem  6  indicates  that  when  P*(k)  is  complete,  k  is  r.n 
upper  bound  for  the  polarized  function  g  on  the  set  P  .  Theorem  8 
naw  states  that  this  upper  bound  may  only  be  attained  on  the  "boundary" 
of  the  polarcid  P*(k)  .  Thus  completeness  means  that  no  interior 

point  of  P*(k)  is  an  optimal  polar  program  x  e  P  ,  with 
k  =  max  g(x)  =  g(x)  .  This  is  stated  in  the  following 


If  P*(k)  is  complete,  then  every  optimal  solution  x  of  the 


irogramming  problem 


maximize  g(xj  ,  subject  tc  .i«:P 
satisfies  x  e  bd  P*(k)  . 

Prop e :  Optimality  implies  x  e  bd  lev^g  >  theorem  8  completes  the  proof. 

Q.E .D. 
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the  above  corollary,  in  a  way,  corresponds  to  the  classical  result 
asserting  that  the  maximum  of  a  quasi-convex  function  over  a  closed 
(bounded)  set  P  is  always  attained  on  the  boundary  of  P  . 

Theorem  9:  Let  Q  c  P  ;  if  P*(k)  is  complete,  then  Q*(k)  is 
complete. 

Proc f :  From  the  inclusion  theorem  1  one  has  P*(k)  ^  Q*(k)  and  com¬ 
pleteness  of  F*(k)  yields 


Q  c  p  <=  levkg  c  (ie\kg)*  (k)  c  p*(k)  c  Q*(k) 

Proof:  Immediate  from  theorem  10.  Q.K.D. 

Theorem  11:  If  (P  vJ  Q)*  (k)  is  complete,  then  both  P*(k)  and 

are  complete. 

Proof:  From  t±  e  union  theorem  one  has 


(P  U  Q)  c  (P  u  Q)*  (k)  =  P*(k)  H  Q*(k)  ; 
but  P  c  (?  !J  Q)  so  that  P  ci  p*(k)  fl  Q*(k)  ; 


and  hence  P  c  P*(k)  ;  and  similarly  for  Q  . 


Q.E.D. 


In  integer  programming ,  one  is  really  interested  in  the  polaroid 


set  of  isolated  sets  (points  in  all  integer-,  or  linear  fibers  in 
mixed  integer -programming ) .  Completeness  of  the  outer-domain  is 
clearly  a  desired  feature  in  order  to  generate  a  deep  cut  and  theorem 
11  indicates  that  such  an  outer  domain  can  be  constructed  as  the  inter¬ 
section  of  individual  polaroids  (one  for  each  point  or  fiber)  proviued 
all  such  polaroids  contain  all  feasible  integer  points. 

The  same  argument  holds  true,  of  course,  in  polar  programming 
if  the  feasible,  set  P  consists  of  (or  is  arbitrarily  split  into) 
several  components. 

5)  Validity  and  optimality  conditions 

The  definition  3  of  a  valid  cut  allows  one  to  consider  optimality 
conditions  in  polar  programming  as  a  particular  valid  cut;  this  is  formu¬ 
lated  in  the  next 

Theorem  12:  The  polarized  function  g(x)  =  f(x,x)  attains  its  maxi¬ 
mal  value  k  =  g(x)  over  the  set  P  at  the  point  xeP 

iff  there  exists  a  valid  cut  S  at  the  level  k  =  max  g(x)  —  g(x) 

xeP 

Proof:  'Trivial  (take  for  instance  S  =  lev_g)  Q.E.D. 

-  k 

In  this  section  we  are  not  only  interested  in  stating  necessary  and 
sufficient  optimality  conditions  for  polar  programs  in  terms  of  more 
general  *and  computationally  more  easily  tractable)  polaroids  than 
the  (trivial)  level  sets  lev  g  ;  we  also  want  to  find  the  conditions 
which  must  be  satisfied,  for  an  arbitrary  cut  to  be  valid;  for  instance, 
the  term  valid  cut  can  here  be  visualize]  as  :  turning  from  the  cutting 
plane  approach  of  Tuy  [7]  for  concave  programming,  or  from  the  inter¬ 
section  cut  approach,  in  integer  programming^  [1,4,5].  In  tne  latter  case,  how¬ 
ever,  a  somewhat  stronger  concept  for  valid  cuts  is  necessary  (see  section  7). 


Theorem  13:  If  the  set  (P  H  S)*  (k)  is  complete  then 
S  is  a  valid  cut  at  the  levtl  k  . 

Proof:  Immediate  from  theorem  7,  applied  to  the  set  (P  H  S)  . 

Theorem  14:  Every  subset  S  of  a  polaroid  P*(k)  is  a  valid  cut. 
Proof:  (Eliminate  the  uninteresting  case  where  S  fl  P  =  0  .)  One  has 
P*(k)  =  {y  |  f(x,y)  <  k  ,  VxeP< 
hence  f(x,y)  <  k  VxeP  ,  VyeS  c:  F*(k)  and  in  particular 
f(x,y)  <  k  ,  Vx,ye  P  fl  S 

Thus  the  function  f  is  complete  on  the  set  (P  fi  S)  ;  theorems  6  and 


7  establish  that  the  definition  3  is  satisfied. 


Q.E.D. 


Corollary  14.1:  P  n  p“(k)  c  levkg 

Proof:  1)  One  may  simply  set  S=P*(k)  and  apply  Theorem  14  together  with 
-  Definition  3. 

2)  Alternately,  the  proof  can  be  obtained  as  follows: 
l‘C(k)  =  fy|f(x,y)  <  k,  VxeP} 

hence  ysP*(k)  f(x,y)  <  k,  VxeP  and,  in  particular, 
when  x=yeP,  one  has 


f  (y,y)  =  g(y)<  k,  VyeP  f.  P  (k) 


Q.E.D. 


In  theorems  13  and  14  one  notes  that  neither  of  than  requites 
completeness  of  ?*(k)  ;  all  that  is  required  in  theorem  13  is  that  the 
cutoff  portion  of  P  (i.e.  (P  fl  S)  )  be  contained  in  the  polaroid 
(p  n  S)*  (k)  ;  from  the  intersection  theorem  3  one  has 


(P  D  S)*  <k)  =>  F*(k)  U  S*  (k) 


I'ITT.TiTV 


showing  that  the  hypothesis  of  theorem  13  is  quite  weak  and  should 
be  easy  to  test  computationally;  in  particular  when  one  constructs  S  as 
a  subset  of  P*(k)  (as  is  often  practically  the  case)  theorem  14 
shows  that  one  merely  must  ensure  S  c  P*  . 

Corollary  6.2:  (Optimality  theorem) 

The  point  x  is  an  optimal  solution  of  the  polar  programming 


problem 


maximize  g(x)  ,  subject  to  xeP 


iff  there  exists  a  complete  polar oid  P*(k)  such  that 


x  e  bd  P*(k) 


Proof:  Immediate  from  theorems  7  and  8. 


Q.E.D. 


6)  Polaroid  cuts  for  linearly  constrained  polar  or  integer  programming  problems . 
T.et  us  now  focus  our  attention  on  the  polyhedral  sets  P  c  i-e. 

P  =  fx  =  x,,...,x  )  I  x.  =  x.  -  j;  a.  .  t.  >  0  ,  Vie (N  U  M)  (2) 

1  n  '  1  1  jeN  13  J  ~ 

where  it  is  assumed  that  (2)  represents  a  linear  program  in  explicit 

(Tucker)  format,  expressed  with  respect  to  a  basis  with  non-basic  set 

N  c  (B  U  H)  ;  N  =  {l.2,...,n}  is  the  set  of  the  original  variables; 

M  =  f n+1 , . . . ,n+m}  Is  the  set  of  the  slack  (and  artificial) variables. 

Consider  the  extreme (basic)  ray  u^(t^)  =  x  -  ajtj  >  t^  >  0  ,  jeN 


where 


Jj  (aij’a2j’**-’anj)  ’ 


x  (xL ,  X,  )  > 

and  assume  (for  simplicity  of  the  exposition)  that  x  e  Int  P*(k)  , 
i.e.  that  the  intersection  points  uJ(\^)  of  each  ray  u~',  jeN  with 
the "boundary  of  P*(k)  are  different  from  x  ,  that  is 

—  7.  rJ:  «-K  o 

U  \A  /  A  *  /v  Uj  }  .Vj.Ln  ..  yj 
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Furthermore  define  the  halfspaces 


H+  =  {x  =  x(t)  j  T,  t  >  1} 

jeN  J 


H*  =  {x  =  x(t)  |  £_  K  t  <  1} 

jeN  3 


where  x(t)  =  (x, (t),x_(t),.. . ,x  (t))  is  given  by  the  n  first  com- 
1  z  n 

ponents  of  the  linear  program  (2)  characterizing  P  .  Clearly  H+  and 
H  are  open,  and  they  are  both  defined  by  the  cutting  plane 

T.  -  t.=l 

jeN  XJ  3 

•  -  -  -f 

(Note  that,  (by  assumption)  xeH  ,i.e.  x  |!  H  ) 

Define  the  set  S  C  P*(k)  as  the  following  n-dimensional  simplicial 

hull  S  =  conv  {x,  uj(Xj)  ,  Vjefl] 

Assuming  that  the  hypothesis  of  one  of  the  theorem  13  or  14  is 
satisfied  the  set  S  represents  a  valid  cut.  (See  section  7)  Implementation  into 
the  current  linear  programming  tableau  merely  amounts  to  the  addi¬ 
tion  of  the  new  linear  constraint 

r_  7.  t.  >  i 
jeN  XJ  3 


called  Polaroid  cut. 


For  the  use  of  a  cutting  plane  method  to  solve  non-convex  prob* 


lems  the  reader  is  referred  to  the  papers  by  Hoang  Tuy  [ 7  ] ,  Glover  [5] 
and  Klingman  f  11  ]  ,  Gomory  [  6  ]  >  Balas  [1,2]  or  Burdet  [  3 ,13].  A  more 


detailed  study  of  particular  types  of  polaroid  cuts  in  integer  pro¬ 
gramming  is  given  in  Burdet  [  4 , 12 ] . 
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7)  Conclusions 

In  polar  programming  the  objective  function  is  the  polarized 
function  g  ;  algorithms  for  solving  such  non-convex  optimization 
problems  therefore  follow  a  stepwise  construction  of  a  set  S 
satisfying 

P  c  s  c  p*(k) 

where  k  is  always  the  current  best  value  and  ultimately 

k  =  max  g(x)  .  The  optimality  theorem  provides  for  vhe 
xeP 

stopping  criterion.  In  practice,  the  use  of  polaroids  will  prove 
efficient  whenever  the  set  P*(k)  is  much  larger  than  P  ;  in  this  case 
ample  room  is  left  for  an  easy  construction  of  a  set  S  yielding  a  suf¬ 
ficiency  condition  for  global  optimality. 

In  integer  programming,  the  use  of  polaroids  is  somewhat  dif¬ 
ferent  since  they  are  only  used  for  the  characterization  of 
the  set  of  feasible  integer  solutions,  as  compared  to  the  other 
(continuous)  feasible  solutions;  the  polarized  function  g  then 
usually  has  the  property: 

g(v)  =1  Vv  =  vertex  of  a  unit  cube  U(x) 

.ontaining  the  linear  programming 
optimum  x  . 

g(x)  <1  Vx  e  U(x) 

thus  k  =  1  ;  the  meaning  of  a  valid  cut  S  here  becomes:  Int(S)  con¬ 
tains  no  feasible  integer  point  (as  shown  by  theorem  13,  the  con¬ 


struction  of  a  valid  cut  requires  neither  completeness  nor  convex- 


"  ‘^-■W i%pVcyizs??&* 
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ity  of  the  Polaroid  P  ).  However,  a  stronger  definition  of  valid  cuts 
is  needed  here;  we  need  add  the  requirement 

Int(S)  n  fp  n  bd  lev^g]  =  0 
to  the  definition  3. 

In  terms  of  this  new  definition.  Theorem  13  remains  true  if  one  requires 
the  additional  hypothesis: 

•k 

SC  (p  n  S)  (k)  ;  indeed,  by  uneorem  8: 

(p  n  s  n  bd  ievkg)  c  (p  n  s  n  bd  (p  n  s)*(k)  > 

every  feasible  integer  point  which  lies  in  S  also  lies  on  the  boundary  of  S  : 

bd  S  3  (p  n  S  H  bd  levkg) 

* 

(Proof  :  since  Int(S)  n  bd  (P  f.  S)  (k)  =  0  one  has 
Int(S)  n  bd  levkg  C  0  ,  Q.E.D.) 

This  general  result  opens  a  new  area  of  research  which  contrasts  from  the 

algebraic  approach  of  Gomory  [6]  or  the  intersection  approach  of  Salas  [1] 

where  one  imposes  in  advance  conditions  which  are  sufficient  to  ensure  the 

validity  of  the  cut  independently  of  the  particular  simplex  S  which 

is  actually  generated.  The  Polaroid  cut  approach  also  entails  such 

a  possibility  (using  convex  and  complete  polaroids) ;  but,  in  addition, 

★ 

it  can  lead  to  new  methods  where  an  arbitrary  polaroid  P  is  used  to 

★ 

generate  the  simplex  S  (owing  to  the  generality  of  P  ,  this  cut 
can  be  made  deep);  in  a  second  phase,  the  validity  of  S  is  then 
established;  for  instance  by  checking  directly  SC  p  (theorem  14). 
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